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ABSTRACT 

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been 
developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and 
mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress 
dependent multifactor interaction relationship developed at NASA Lewis Research Center has been used to model 
the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to 
compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive 
random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed 
and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the 
thermal cyclic loads on the fatigue reliability for a (0/±45/90) s graphite/epoxy laminate with a ply thickness of 
0.1 27 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical 
cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix com- 
pressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical 
cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear 
strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness. 


INTRODUCTION 

Under the NASA High Speed Research (HSR) Program, technological solutions which will ensure the 
economic viability and environmental compatibility of a future High Speed Civil Transport (HSCT) plane are cur- 
rently being sought. Lighter structural concepts for both airframe primary structures and engine structure com- 
ponents are being investigated. It is envisioned that such objectives can be achieved through the use of high 
temperature composites as well as other conventional lighter weight alloys. One of the prime issues for these struc- 
tural components is assured long term behavior with a specified reliability. For example, the engine structure com- 
ponents are required to last 18 000 hr, and the airframe structural components are required to last 60 000 hr of flight 
time. These objectives are being pursued with the cooperation of industry and university participants through 
developments in accelerated testing strategies as well as computational simulation methodologies for predicting the 
long term durability of HSCT structural components. The present report addresses issues pertaining to probabilistic 
fatigue life of polymer matrix composites under thermal/mechanical cyclic loads through some typical examples. 

The complexity of predicting composite behavior is compounded by the presence of multiple scales that 
must be considered in the simulation. These are the micro level scale where fiber, matrix, and interface participate, 
the macro level where individual plies take part, and the laminate level where the ply configuration becomes impor- 
tant. Additionally, the fabrication process induced variations, inherent uncertainties in the constituent properties, and 
cyclic variations in an aggressive loading environment all contribute to the ever increasing complexity in the compu- 
tational simulation. Successful inclusion of the lighter high temperature composites in future aircraft hinges largely 
upon the ability to predict and assure their long term behavior. Current practices depend largely on a building block 
approach of conducting long term testing of materials in the form of coupons to components. However, it is time 
consuming and expensive to conduct a large number of long term tests to capture the effects of all the design 
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variables and loads. Furthermore, accounting for the previously discussed uncertainties in the experiments is mon- 
etarily unfeasible proposition. Therefore in view of these facts, realization of the full potential of lightweight com- 
posites is more likely to be achieved with the real innovative computational approaches capable of handling these 
aspects in an integrated manner together with some accelerated and strategically chosen experiments to minimize the 
testing costs. Also such computational procedures provide a priori insight to identify critical experiments to be per- 
formed as well as reduce the number of tests to be conducted to an optimum number. 

Traditional computational approaches are deterministic in nature; they do not account for uncertainties 
associated with composite structures and materials. The focus of ongoing research at the NASA Lewis Research 
Center has been to develop advanced integrated computational methods and related computer codes to perform a 
complete probabilistic assessment of composite structures. These methods account for uncertainties in all the con- 
stituent properties, fabrication process variables, and loads to predict probabilistic micro, ply, laminate, and struc- 
tural responses. These methods have already been implemented in the Integrated Probabilistic Assessment of 
Composite Structures (JPACS) (ref. 1) computer code. 

The objective of this report is to describe a computational simulation methodology to predict the fatigue 
life reliability and probabilistic long term behavior in polymer matrix composites. Also, the implementation of the 
methodology in the ICAN computer code (ref. 2) has been described. A unified time, stress, and load dependent 
multifactor interaction equation (MFIE) model (ref. 3) developed at the NASA Lewis Research Center has been 
used to simulate the long term behavior of polymer matrix composites. 

To illustrate the application of the methodology, a typical composite system consisting of graphite fibers 
in epoxy matrix is chosen. The methodology can be applied to other types of polymer matrix composites as well. 

The layup of the composite system is chosen as (0/±45/90) s . The cumulative probability distribution functions 
(CDF's) for the fatigue life cycles are computed for different applied stress to laminate strength ratios, r. The lami- 
nate strength is evaluated based on first ply failure criteria (hereinafter referred to as laminate strength). First ply 
failure criteria assumes that the laminate has failed when any stress component in a ply exceeds its respective allow- 
able limit. By using these CDF's a fatigue life cycle curve for a reliability of 0.999 is obtained in order to demon- 
strate how this methodology can be used to aid the design process. Also, the sensitivity of fatigue life reliability to 
the primitive random variables are computed for a reliability of 0.999. 
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NOMENCLATURE 

fiber normal modulus 

matrix normal modulus 

fiber shear modulus 

actual fiber volume ratio 

apparent void volume ratio 

exponents for multifactor interaction equation 

material property 

reference material property 

number of effects causing material degradation 

number of mechanical cycles 

final number of mechanical cycles 

number of thermal cycles 
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N Tf final number of thermal cycles 

p coefficient of variation 

r applied stress to laminate ratio 

S final strength 

Sf fiber strength 

Sfc fiber compressive strength 

Sfr fiber tensile strength 

S m c matrix compressive strength 

S m s matrix shear strength 

S m T matrix tensile strength 

T temperature 

T g d dry glass transition temperature 

Tgw wet gears transition temperature 

T 0 initial temperature 

t current time 

tf final time 

ti ply thickness 

V primitive variable for material or load (e.g., temperature, stress, mechanical cycles, etc.) 

Vp value of primitive variable at final condition 

Vo value of primitive variable at reference condition 

0 Cm matrix thermal expansion coefficient 

ji mean value 

Vfi 2 , 13,23 fiber Poisson’s ratio 

v m i2, 13,23 matrix Poisson’s ratio 

o stress 

Cm mechanical cyclic stress amplitude 

Gt thermal cyclic stress amplitude 
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COMPUTATIONAL SIMULATION 


Simulation of uncertainties in the primitive random variables throughout the load history requires 
fundamental governing field equations in time variable. The probabilistic simulation process in this report uses 
time dependent MF1E model, micromechanics theory, and the fast probability integration (FPI) theory (ref. 4). 
Uncertainties in the primitive variables related to physical process fabrication, constituent properties, and ply geom- 
etry are included in the simulation. Micromechanics theory available in the ICAN computer code has been used 
together with the MFIE to include the degradation of constituent properties resulting from aggressive load effects for 
the simulation of long term effects. A generic form of the previous model is given by 
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The terms in parenthesis represent a physical effect that causes the degradation of material property. The exponent 
accounts for the amount of degradation due to a particular effect. More details on this model are given later in this 
report. The following sections describe the developed/modified computer codes and procedures used in the com- 
putational simulation of probabilistic long term behavior in polymer matrix composites. 


Integrated Composite Analyzer (ICAN) 

The ICAN code computationally simulates the material behavior of polymer matrix composites from fiber/ 
matrix constituents to the laminate scale including fabrication effects (fig, 1). ICAN uses composite micromechanics 
and laminate analysis based on linear elastic theory to compute constituent, ply, and laminate level properties 
required for the global structural analysis (left side of fig. 1). ICAN also decomposes the global structural response 
to laminate, ply , and constituent response levels (stress and strains) which helps the user to evaluate failure (right 
side of fig. 1). Details of the ICAN computer code are given in reference 2. ICAN also performs the failure analysis 
based on different failure criteria such as first ply failure criteria, fiber break criteria, modified distortion energy 
criterion, and Hoffman’s criteria. 

The ICAN computer code was modified to implement the time dependent MHE and perform sensitivity 
evaluation for primitive random variables. The time variable in the MFIE is divided into a smaller number of time 
steps to simulate long term effects. At every time step the degradation magnitude of the matrix properties, which are 
used in micro-mechanics and laminate analysis, is evaluated. Response variable sensitivity to the random variables 
at every time step are also performed to compute the respective scatter in response variables. Also, at every time step 
a failure analysis based on first ply failure and fiber break criteria is performed to determine whether the laminate 
can take any further load. Failure analysis determines the possible failure modes and maximum load capacity in the 
respective failure mode. The analysis stops when the laminate is incapable of carrying any more load. In addition, a 
data base is created to store results required to perform probabilistic life analysis and develop reliability based life 
assessment of polymer matrix composites. Thus, the current version of ICAN is capable of performing probabilistic 
time variable analysis to incorporate material degradation due to thermal and mechanical cyclic loads. 


Time Dependent Multi-Factor Interaction Equation (MFIE) Model 

It is known that predicting the behavior of composite materials is a complex task. Accounting for all the 
physical effects and how they affect the material properties in time domain is even more complex. Over the years, 
research in developing a unified law describing the material behavior driven by primitive variables has been an 
ongoing activity at NASA Lewis Research Center. The result of this research is the development of a unified MFIE 
model (ref. 3). Concepts used in reference 3 have been expanded to include the time dependent degradation effect on 
material behavior because of environmental, fabrication, and load effects (ref. 5). The modified form of the equation 
including time effect is 
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Each term in parenthesis accounts for a specific physical effect. Any number of effects can be included in one single 
equation as seen by the nature of the equation. The exponents are determined from the available experimental data 
or estimated from the anticipated material behavior due to a particular primitive variable. Each primitive variable 
and the exponent in the previous equation can be random with a statistical distribution. The insufficiency of a set of 
experimental data can be taken into account by means of uncertainties in the exponent. 

An important part of the previous model is the fact that only one equation can include all the effects with 
any nonlinearity in the material behavior and follow the physics of behavior. It can describe all the interacting ef- 
fects of different variables (thermal, metallurgical, mechanical, chemical, and load). Since the variables used are at a 
primitive level, it simulates the in situ degradation in material properties caused by applied cyclic and environmental 
effects. 


Probabilistic Simulation 

An advanced first-order second-moment fast probability integration (FPI) technique (ref. 4) is used to com- 
pute cumulative probability distribution function of the fatigue life. Results of the random variable perturbations are 
used to compute the CDF of responses. The FPI technique is very efficient compared to the Monte Carlo simulation 
technique. FPI does not generate random samples; however, it does use a numerical integration technique to com- 
pute the joint probability density functions and probability of failures. It transforms the physical random variable 
space system into unit normal space to perform probability integration easily and more accurately. On the 
other hand, the Monte Carlo method generates a large number of random samples to compute the CDF of response. 
Hence, FPI is computationally more economical than Monte Carlo. The sensitivity of responses is also computed by 
FPL Sensitivity information helps improve the design and quality. 


SIMULATION CASES, RESULTS AND DISCUSSION 

Probabilistic simulation of long term behavior under mechanical cyclic loads only was demonstrated in 
reference 5. In this report the methodology has been extended to simulate probabilistic long term behavior and 
fatigue reliability under combined thermo-mechanical cyclic loads. Probabilistic behavior of a graphite/epoxy 
laminate as an illustrative example is studied. The deterministic long term behavioral aspects for polymer matrix 
composites are discussed in reference 6. A (0/±45/90) s laminate made of graphite fibers and epoxy matrix is sub- 
jected to uniaxial mechanical tensile and thermal cyclic load as shown in figure 2. The load shown in figure 2 is 
sinusoidal. However, any cyclic load shape such as triangular, square, and so forth can be used by accounting for it 
through the exponent value in MFIE. Each ply has a thickness of 0.127 mm. Typical mean values, scatter, and prob- 
ability distribution of all the fiber, matrix properties, and fabrication variables are taken as those given in table I (see 
ref. 2). Normal distribution for the primitive variables has been assumed in the present study. However, any other 
type of distribution for these variables can be used in the described methodology. Initially, the laminate was sub- 
jected to a mean static load, and failure analysis was performed to evaluate static laminate strength. It was found to 
be 1184.0 kg per unit cm width based on the first ply failure criteria. Several simulation cases with different combi- 
nations of mechanical and thermal load magnitudes given in table II were analyzed, and the CDFs and sensitivity 
evaluations were performed. The mechanical load is applied as a ratio, r, defined as applied mechanical stress over 
static laminate strength. These load cases are divided into three categories. The first category includes only the 
mechanical cyclic load (load cases 1 to 3 in table II). However, the category 2 load cases are only for thermal cyclic 
loads (load cases 4 to 7 in table 13). In the category 3 load cases, the mechanical cyclic loads are combined with the 
thermal cyclic loads of different amplitudes such as 0, 15.6, 23.9, 37.8, and 51.7 °C (load cases 8 to 22 in table II). 
The purpose of having these categories is to illustrate the combined effect of thermal and mechanical cyclic loads on 
the fatigue life reliability. 
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In each load case simulation, deterministic sensitivity evaluation of the fatigue life to each random variable 
was performed using perturbation analysis (ref. 7). Perturbation analysis means evaluating the change in the 
responses if the value of a particular random variable is changed by a given magnitude of standard deviation of that 
variable. From each perturbation analysis, the fatigue life for each ply was computed based upon the first ply failure 
criteria. The shortest life of a ply in the laminate was considered to be the life of a laminate for that perturbation. 
Perturbation analysis is performed to develop a response surface which is required for the probabilistic analysis. 

By using these perturbation analysis results and probability distributions of random variables, the reliability analysis 
was performed using the FPI method. The CDF’s and sensitivity factors for the first three load cases are summarized 
in figure 3. It shows the CDF of fatigue life caused by mechanical cyclic loads only. The fatigue life in these three 
cases was primarily governed by failure in 90° ply because the transverse tensile stress was greater than the cor- 
responding strength. Thus, the life under uniaxial tensile cyclic loads of that laminate is governed by the transverse 
tension failure in 90° plies. Figure 3 also shows that the fatigue life CDF shift to the left means that the life under 
mechanical cyclic loads drops rapidly as the amplitude of applied cyclic stress increases. At load ratios, r, equal to 
60, 70, and 80 percent, the respective mean fatigue life, ji, was computed to be 88, 40, and 16 percent of the “endur- 
ance limit, N^.” The corresponding scatter, p, with respect to mean value is 14.4, 26, and 32 percent. Also, the 
respective fatigue life for 0.999 reliability is 49, 9.3, and 0.6 percent of as shown in figure 4 (S-N curve) and 
table II. Thus, the fatigue life drops rapidly and the p increases. Since ply transverse tension is the governing failure 
mode, the fatigue life for 0.999 reliability is sensitive to the S mT , ply thickness, E m , E fl v and ky. Note that the 
magnitude of sensitivity for a particular random variable almost remains the same with respect to different values of 
r. This means that the scatter in these variables should be minimized or controlled in order to enhance the fatigue life 
reliability under mechanical loads. 

Figure 5 depicts the CDF’s and random variable sensitivity for the fatigue life caused by variation in the 
cyclic temperature amplitude (load cases 4 to 7 and table II). Fatigue life caused by thermal cyclic loads in the lami- 
nate under consideration is mainly governed by the transverse compression failure in 90° and 45° plies. It is seen 
from this figure that the fatigue life variation is not as severe due to variation in cyclic temperature amplitude as was 
the case for mechanical load amplitude. The computed mean fatigue life, m, for temperature amplitudes of 15.6, 
23.9, 37.8, and 51.7 °C is 91, 86, 79, and 73 percent of N^, respectively, whereas the scatter, p, is 6.47, 6.91, 6.57, 
and 6.73 percent, respectively. Also, the respective fatigue life for 0.999 reliability is 73, 68, 63, and 58 percent 
(fig. 6 and table El). Thermal loads induce compressive stresses in the transverse direction of both 45° and 90° plies. 
Since allowable ply compressive strength is higher than the tensile strength, a higher life is expected under thermal 
loads. Also, as seen from figure 5, E m , a m , S mC , k^ S mS , T gd , ply thickness, and ky govern the fatigue life for 0.999 
reliability. Again, the magnitude of sensitivity for different thermal cyclic load amplitudes does not change much. 
The sensitivity trend of fatigue life to S mC is downward with the increase in the thermal load amplitude which is 
indicative of other eminent failure modes at very high thermal amplitudes. It is interesting to note that the scatter, p, 
remains the same at all amplitude levels. 

Figures 7 to 1 1 are CDF’s and sensitivity factors for the fatigue life under combined thermo-mechanical 
loading with cyclic temperature amplitudes of 0, 15.6, 33.9, 37.8, and 51.7 °C, respectively. Figure 12 shows the 
fatigue life design curves for 0.999 reliability. The mean fatigue life, scatter, and the life for 0.999 reliability for these 
load cases are given in table H. It is seen from figures 7 and 8 that the fatigue life CDF shifts to the left in proportion to 
the applied load ratio, r, at low thermal amplitudes of up to 23.8 °C which means that the probability distribution of 
fatigue life does not change much. This is due to the fact that the failure mode up to 23.9 °C cyclic amplitude remains the 
same. This results in linear variation of the mean, scatter, and 0.999 reliability value of fatigue life (fig. 12). The vari- 
ables that control the fatigue life magnitude at 0,999 reliability are ply thickness, S mC , S mS , E m , a m , T gd , kf, ky, and 
S mT . Ply thickness is the important variable in these cases caused by the mechanical load effects. 

As the cyclic thermal amplitude increases, the shift in fatigue life CDF’s is not in proportion to r as seen 
from figures 9 to 1 1 . Thus, the trend in variation of m, p, and 0.999 reliability value of fatigue life gradually be- 
comes nonlinear implying that the nature of probability distribution of fatigue life is not same at all load ratios. For 
amplitudes of 37.8 and 51.7 °C it becomes highly nonlinear (refer to fig. 12 and table II). This is due to the fact that, 
at low temperature amplitudes, the failure is mostly governed by the longitudinal compressive strength in a 90° ply 
whereas, at higher thermal amplitudes and low load ratio, the governing failure mode shifts from 90° ply to 45° ply 
with compressive strength in the transverse direction being the limiting one. Also, at high load ratio and high tem- 
perature amplitudes, the failure mode is controlled by the shear strength in a 45° ply. Thus, when the shear mode 
failure controls, the drop in fatigue life is also steep and the scatter in fatigue life increases dramatically. As 
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evidenced by sensitivity factors depicted in figures 9 to 11, the fatigue life at 0.999 reliability at low load ratio is 
dominated by S mC , E m , a m , kf, and S mS . However, at higher load ratios and thermal amplitudes, the fatigue life sen- 
sitivity to S mS , E fll , and ply thickness increases whereas the others’ effect decreases. Thus, it is seen that the nature 
of failure mode governs the reliability and corresponding variables become sensitive to the life. Also, it is observed 
that the shear failure mode causes a large scatter in the fatigue life as well as reduction in the life for 0.999 reliabil- 
ity. It can be seen from figure 12 that for the laminate under consideration, the 0.999 reliability cannot be achieved 
for load ratio, r, of 0.9 and thermal amplitudes of 37.8 and 51.7 °C. At this magnitude of reliability the fatigue life is 
negative in figure 12, meaning the laminate has already failed at this reliability. 

The information in figures 4, 6, and 12 can be used as a design aid to determine the component life for 

0.999 reliability as well as to make decisions on its continuing use. Similar design curves can be developed readily 
for any polymer matrix composite with any laminate configuration and for any load conditions. It is also important 
to mention here that, when the failure mode changes for different perturbations under same load conditions, one 
should iterate the probabilistic analysis to ensure the convergence of reliability using advanced first order mean 
value analysis given in reference 7. This is necessary because a small perturbation in the random variables could 
change the entire failure mechanism and the response surface becomes highly non-linear which in turn affects the 
probability of failure. A reliability-based approach provides the level of reliability of a given design as well as iden- 
tifies the controlling variables. Thus, it helps improve the manufacturing process, design experiments, and develop 
inspection guidelines. 


CONCLUSION 

A methodology to simulate the probabilistic fatigue life of polymer matrix composites under thermal and 
mechanical loads has been developed, implemented in the in-house computer code Integrated Composite Analyzer 
(ICAN), and demonstrated by the examples. The methodology incorporates the ICAN computer code, a generic 
time dependent multi-factor interaction equation model, and fast probability integration technique. Cumulative dis- 
tribution functions of fatigue life cycles for a (0/±45/90) s graphite/epoxy laminate under uniaxial tensile cyclic load 
and thermal cyclic loads were computed. Sensitivities of fatigue life to the random variables for 0.999 reliability 
were plotted. Fatigue life design curves for 0.999 reliability were generated for different combinations of thermal 
cyclic load amplitude and mechanical cyclic load. The importance of the failure mode on fatigue life reliability, 
scatter, and sensitivity were discussed. The fatigue life for (0/±45/90) s graphite/epoxy laminate under mechanical 
cyclic loads only was most sensitive to thickness and matrix tensile strength at 0.999 reliability since tension failure 
in 90° ply controls the life. However, at low mechanical cyclic loads and low thermal cyclic amplitudes, the matrix 
compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness were most sensitive to the 
fatigue life for 0.999 reliability since the transverse compression failure in 45° and 90° plies controls the life. Also, at 
high mechanical cyclic loads and high thermal cyclic amplitudes, shear failure in 45° ply controlled the life. There- 
fore, at these load magnitudes, the fatigue life at 0.999 reliability was more sensitive to the shear strength of matrix, 
longitudinal fiber modulus, matrix modulus, and ply thickness. The reliability-based fatigue curves are useful in 
determining the allowable fatigue load on the structure or assessing the life of a given component. Sensitivity infor- 
mation provides guidelines to the designer to improve the reliability of the structure’ s fatigue life. 


REFERENCES 

1. Integrated Probabilistic Assessment of Composite Structures (IPACS) User’s Manual. Ver. 1.0, NASA Contract 

No. NAS3-25266, Task Order No. 5209-4)2, Dec. 1991. 

2. Murthy, P.L.N.; and Chamis, C.C.: Integrated Composite Analyzer (ICAN): Users and Programmers Manual, 

NASA TP-2515, March 1986. 

3. Boyce, L.; and Chamis, C.C.: Probabilistic Constitutive Relationships for Cyclic Material Strength Models. 

Proceedings of the 29th ALAA Structures, Structural Dynamics, and Materials Conference, pt. 3, J. Propul. P., 
1988, pp. 1299-1306. 

4. Wu, Y.T.: Demonstration of a New Fast Probability Integration Method for Reliability Analysis. Advances in 

Aerospace Structural Analysis, O.H. Burnside and C.H. Parr, ed., ASME, New York, pp. 63-73, 1985. 


7 



5. Shah, A.R.; Singhal, S.N.; Murthy, P.L.N.; and Chamis, C.C.: Probabilistic Simulation of Long Term Behavior in 

Polymer Matrix Composites. Proceedings of the 35th AIAA/ASME/AHS/ASCE Structures, Structural 
Dynamics, and Materials Conference, April 1994. 

6. Shah, A.R; and Chamis, C.C.: Cyclic Load Effects on Long Term Behavior of Polymer Matrix Composites. 

Proceedings of the 40th International SAMPE Symposium and Exhibition, May 1995. 

7. Probabilistic Structural Analysis Methods (PS AM) for Select Space Propulsion System Components User’s 

Manual, South West Research Institute, NASA Contract NAS3-29389, Nov. 1989. 


TABLE L— UNCERTAINTIES IN THE 
CONSTITUENT PROPERTIES 
[Scatter, 5 percent; normal distribution.] 


Property 

Mean 

Fiber 

Normal modulus, E fll , GPa 

213,7 

Normal modulus, E G2 , GPa 

13.8 

Poisson's ratio, v fl2 

G.20 

Poisson’s ratio, 

0.25 

Shear modulus, G fl2 , GPa 

13.8 

Shear modulus, G^, GPa 

6.9 

Tensile strength, S^, GPa 

2.8 

Compressive strength, S^, GPa 

2.8 

Matrix 

Normal modulus, E m , GPa 

3.4 

Poisson’s ratio, v ml2 

0.35 

Tensile strength, S mX , GPa 

0.1 

Compressive strength, S mC , GPa 

0.24 

Shear strength, S mS , GPa 

0.09 

Fabrication variables 

Fiber volume ratio, k f , percent 

60 

Void volume ratio, k v percent 

2 

Ply thickness, tj , mm 

0.127 
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TABLE H.— COMPUTATIONAL LOAD CASES AND 
RELIABILITY-BASED FATIGUE LIFE 


Case 

Mechani- 
cal cyclic 
load ratio, 
r 

Mean 
temperatu 
re, °C 

Cyclic 

temperature 

amplitude, 

°C 

Mean 
fatigue life 

XNmf 

Scatter, 

percent 

0.999 

Reliability 

life 

xN tf 

1 

0.6 

21.1 

0.0 

0.8792 

14.40 

0.4909 

2 

0.7 

21.1 

0.0 

0.3956 

25.65 

0.0934 

3 

0.8 

21.1 

0.0 

0.1609 

32.3 

0.0068 

4 

0.0 

65.6 

15.6 

0.9080 

6.47 

0.7318 

5 

0.0 

65.6 

23.9 

0.8623 

6.91 

0.6788 

6 

0.0 

65.6 

37.8 

0.7909 

6.57 

0.6287 

7 

0.0 

65.6 

51.7 

0.7254 

6.73 

0.5758 

8 

0.6 

65.6 

0.0 

0.9045 

5.1 

0.7616 

9 

0.7 

65.6 

0.0 

0.7785 

5.42 

0.6486 

10 

0.8 

65.6 

0.0 

0.6423 

6.5 

0.5118 

11 

0.6 

65.6 

15.6 

0.9278 

5.6 

0.7529 

12 

0.7 

65.6 

15.6 

0.8076 

5.26 

0.6744 

13 

0.8 

65.6 

15.6 

0.6736 

4.03 

0.5516 

14 

0.6 

65.6 

23.9 

0.8980 

6.63 

0.7126 

15 

0.7 

65.6 

23.9 

0.8231 

5.51 

0.6764 

16 

0.8 

65.6 

23.9 

0.6898 

6.57 

0.5433 

17 

0.7 

65.6 

37.8 

0.8186 

5.45 

0.6610 

18 

0.8 

65.6 

37.8 

0.6949 

6.0 

0.5590 

19 

0.9 

65.6 

37.8 

0.4009 

44.0 

-0.1154 

20 

0.7 

65.6 

51.7 

0.7709 

6.72 

0.6143 

21 

0.8 

65.6 

51.7 

0.7047 

6.15 

0.5710 

22 

0.9 

65.6 

51.7 

0.3984 

44.0 

-0.1159 
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Figure 1. — Schematic of ICAN computer code. 
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Figure 2. — Description of thermo-mechanical cyclic load. 
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Figure 3.— Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite 
epoxy laminate. Ply thickness, 0.1 27 mm under mechanical cyclic load; mean mechanical load = 0.5 x static strength, 
(a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 4.— Fatigue life variation for 0.999 reliability of 
(0/±45/90) s graphite/epoxy laminate. Ply thickness, 
0.127 mm under mechanical cyclic load. Mean 
mechanical load = 0.5 x static strength. 
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Figure 5. — Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite/epoxy 
laminate. Ply thickness, 0.127 mm under thermal cyclic load; mean mechanical load 0.5 x static strength, (a) Cumulative 
distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 6.-— Fatigue life variation for 0.999 reliability of 
(0/±45/90) s graphite/epoxy laminate. Ply thickness, 
0.127 mm under thermal cyclic load; mean mechanical 
load = 0.5 x static strength. 
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Figure 7.— Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite/ 
epoxy laminate. Ply thickness, 0.1 27 mm under thermo-mechanical cyclic load; mean temperature, 65.6 °C; amplitude 
temperature, 0 °C. (a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 8.— Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite/ 
epoxy laminate. Ply thickness, 0.127 mm under thermo-mechanical cyclic load; mean temperature, 65.6 °C; amplitude 
temperature, 1 5.6 °C. (a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 9. — Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite/ 
epoxy laminate. Ply thickness, 0.1 27 mm under thermo-mechanical cyclic load; mean temperature, 65.6 °C; amplitude 
temperature, 23.9 °C. (a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 10. — Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90) s graphite/ 
epoxy laminate. Ply thickness, 0.127 mm under thermo-mechanical cyclic load; mean temperature, 65.6 °C; amplitude 
temperature, 37.8 °C. (a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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(a) N m /N mf (b) 

Figure 1 1 . — Cumulative distribution function for and sensitivity of fatigue life for 0.999 reliability of (0/±45/90)s graphite/ 
epoxy laminate. Ply thickness, 0.1 27 mm under thermo-mechanical cyclic load; mean temperature, 65.6 °C; amplitude 
temperature, 51 .7 °C. (a) Cumulative distribution function, (b) Sensitivity of fatigue life for 0.999 reliability. 
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Figure 12. — Fatigue life variation for 0.999 reliability of 
(0/±45/90) s graphite/epoxy laminate. Ply thickness 

0.1 27 mm under thermo-mechanical cyclic load; mean 
temperature, 65.6 °C. 
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